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(D ■ ABSTRACT 

cn 

I We describe results from a new ground-based monitoring campaign on 
NGC 5548, the best studied reverberation-mapped AGN. We find that it was 

I in the lowest luminosity state yet recorded during a monitoring program, namely 

^ I -Z^5ioo = 4.7 X 10'*^ ergs s~^. We determine a rest-frame time lag between flux vari- 

^ I ations in the continuum and the H/? line of 6.3I2 3 days. Combining our mea- 

O ' surements with those of previous campaigns, we determine a weighted black hole 
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mass of Mbh = 6.541q25 x IO^Mq based on all broad emission lines with suitable 
variability data. We confirm the previously-discovered virial relationship between 
the time lag of emission lines relative to the continuum and the width of the emis- 
sion lines in NGC 5548, which is the expected signature of a gravity-dominated 
broad-line region. Using this lowest luminosity state, we extend the range of the 
relationship between the luminosity and the time lag in NGC 5548 and measure 
a slope that is consistent with a = 0.5, the naive expectation for the broad line 
region for an assumed f orm of r oc L ". This value is also consistent with the slope 



recently determined by iBentz et al.l for the population of reverberation-mapped 
AGNs as a whole. 

Subject headings: galaxies:active — galaxies: nuclei — galaxies: Seyfert 



INTRODUCTION 



Reverberation mapping (iBlandford fc McKedll982l : |Petersorull993l ) is an extremely use- 
ful technique that has been exploited to measure the size of the broad-line region (BLR) 
for many nearby active galactic nuclei (AGNs). The time delay, or lag, between variations 
in the continuum flux and broad line flux (usually H/3) provides a measure of the size of 
the region from which the broad lines emanate. Coupling the measured lag time with the 
velocity width of the emission line provides an estimate of the mass of the central black hole. 
To date, 36 nearby Seyfert g alaxies have BLR radius a nd black hole mass determinations 
from reverberation mapping (jPeterson et al.l |2004J . |2005| ) . 



The Seyfert galaxy NGC 5548 was the focus of an inte nse, 13-year campaign by the 
International AGN Watch consortiunj^ ( jPeterson et al. 2002 ) to study the variations in the 
optical continuum and H/3 line flux. As a result, its emission line variability properties are 
well characterized. 

During the spring of 2005, we un dertook a new ground-based monitoring program 
(IBentz et al.l l2006bl : iDenney et al.l |2006| ) with the principal aim of replacing H/5 lag mea- 
surements and black hole mass estimates for several AGNs with previous, unsatisfactory 
data. We also monitored NGC 5548 and found that it was in the lowest luminosity state yet 
recorded during a monitoring campaign. In this paper, we present the low-luminosity moni- 
toring data and combine it with previous monitoring results of NGC 5548 to examine several 
questions related to reverberation-mapping and AGN variability. We explore the behavior 



http://www.astronomy.ohio-state.edu/~agnwatch/ 
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of the BLR with respect to varying luminosity states in NGC 5548 as well as examine the 
intrinsic scatter in reverberation-based black hole mass measurements for a single object. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Spectroscopy 

For most of the clear nights between 2005 March 1 and 2005 April 10, we obtained 
spectra of the nucleus of NGC 5548 — an SO/a Seyfert galaxy at z = 0.01717 — with the 
Boiler and Chivens CCD Spectrograph (CCDS) on the McGraw-Hill 1.3-m telesc ope at MDM 



Observ atory. The observational details and data reduction are as described by iBentz et al. 



(l2006bl ). but for completeness, we include a few of the details here. The observations of 
NGC 5548 were obtained at a position angle of 90° through a 5" slit, resulting in a spectral 
resolution of 7.6 A over the wavelength range 4400 — 5650 A. The typical seeing was 2". 

Additional spectra were obtained at the Crimean Astrophysical Observatory (CrAO) 
2.6-m Shain Telesco pe with the Nasmith Spectrograph and Astro-550 580 x 520 pixel CCD 



(IBerezin et al.lll99ll ). The spectra were obtained through a 3" slit at a position angle of 90°. 
Spectral reduction was carried out in the usual way with an extraction width of 16 pixels, 
corresponding to 11" on the sky. 

The spectra were then internally flux calibrated within each data set. We scaled each 
individual spectrum to an [O III] A5007 flux of 5.58 x 10^^'^ ergs s~^ cm~^ — the flux of 
the [O III] A5007 line determined by photometric spectra from the fir st year of the 13-year 



ground-based monit oring program for NGC 5548 (IPeterson et al.lll99ll ) — using the spectral 



scaling algorithm of Ivan Groningen fc Wanders! (119921 ). This process minimizes the residuals 



of the [O III] A5007 line in a difference spectrum by comparing individual spectra to a 
reference spectrum (in this case, the mean spectrum of all the data in the set) and scaling 
appropriately. 

We then created light curves from measurements of the spectra. The mean continuum 
flux was measured between observed-frame wavelengths of 5170-5200 A. The flux of the H/? 
emission line was measured by first setting a linear continuum level between the observed- 
frame windows of 4800-4830 A and 5170-5200 A, and then integrating the flux above the 
best-fit continuum between 4850-5000 A in the observed frame. 
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2.2. Photometry 



Photometry in the ^-band was obtained at the 2.0-m Muhicolor Active Galactic NUclei 



Monitoring; (MAGNUM) telescope at the Haleakala Observatories in Hawaii (iKobayashi et al 



1998bj : lYoshiill2002l : lYoshii et al.ll2003l ) with the m ulticolor imaging; photorn eter (MIP) (IKobayashi et al 
1998al ). Following the procedure explained by ISuganuma et al.l (120041 ). the nuclear flux 
of NGC 5548 was measured relative to a nearby reference star located at {Aa,A6) = 
(— Ofl, — 2'6). After processing the frames in the usual way with IRAF, aperture photometry 
was performed within an aperture radius of 4'.' 15, with sky subtraction between radi i of 5^^5 



&!9. T he r eference star was t hen calibrated using photometric standard stars from iLandolt 



Jl992h and lHunt etaP Jl998l) Vari ability of the reference star was previously shown to be 
negligible by ISuganuma et "aL (2004). Due to technical issues with the instrument, the pho- 
tometric accuracy was very low at the beginning of this monitoring campaign, but improved 
considerably towards the end when the instrument was flxed (note the magnitude of the 
uncertainties listed in Table 1 and the size of the error bars in Fig. 1). 



2.3. Inter-calibration of Light Curves 

As a result of the varying slit geometries, extraction apertures, and seeing conditions 
in each of the spectroscopic and imaging data sets, they must be inter-calibrated with each 
other to account for slit l osses and different am ounts of host-galaxy starlight flux. We follow 



the process described by iBentz et al.l ( l2006bl ) to scale the CrAO and MAGNUM data sets 



to the MDM data set by using a least-squares analysis to identify the flux offsets relative to 
the MDM data set. After correcting for these offsets, the three data sets were merged into a 
continuum light curve and B.(3 light curve, which are tabulated in Table 1. For the analysis 
that follows, we combined observations within a 0.5 day window, which yields the final light 
curves that are depicted in Figure 1. The points with the large error bars are those from the 
MAGNUM telescope while the instrument was experiencing difficulties. 

Table 2 describes the statistical properties of the final continuum and emission-line 
light curves. Column (1) lists the spectral feature of the time series and column (2) gives the 
number of measurements. The mean and median sampling interval between data points are 
given in columns (3) and (4), respectively. The mean fiux and standard deviation of the time 
series are given in column (5). The continuum fiux measurement does not take into account 
the host-galaxy starlight contribution. Using a h igh-resolution Hubbl e Space Telescope image 



of NGC 5548 and the procedures described by iBentz et al.l (l2006al ). we measure the host- 
galaxy contribution through the larger slit geometry in this study (5'.'0 x 12'.'75) to be 
Fgai(5100A) = (5.29 ± 0.49) x lO-^^ ergs s ^ cm^ -2 A-i. The mean fractional error, which 
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is based on closely spaced observations in time, is given in Column (6). Column (7) is the 
excess variance, which is computed as 



if) 



(1) 



where is the variance of the fluxes, 5^ is their mean-square uncertainty, and (/) is the 
mean of the observed fluxes. Lastly, column (8) is the ratio of the maximum to the minimum 
flux in each time series. Comparison of the values of Fvar for the continuum and H/3 fluxes 
with those determined for the previous 13 years of monitoring data show that the amplitude 
of variation in the con tinuum is about half that of the previous lowest continuum- variability 
year (Year 7 [1995] of IPeterson et al.ll2002h and the amplitude of variabil ity in the H/3 li ne 
is the third lowest measured (see also Year 5 [1993] and Year 11 [1999] of lPeterson et al.l ). 



Correcting the mean optical flux for the host-galaxy contribution above, we find that 
NGC 5548 was in the lowest luminosity state yet recorded during a monitoring campaign. 



only Fagn(5100A) = (1.3±0.5) x 10"^^ ergs s"^ cm'^ A'^ This is a full 40% fainter than the 



previous record holder. Year 4 (1992) of the AGN Watch 



for the original flux measurement and lBentz et al.ll2006al for the host galaxy correction). 



campaign (see 



Peterson et al. 2002 



3. DATA ANALYSIS 



3.1. Time Series 



As a result of the low luminosity of the AGN and the very low level of variability 
throughout the campaign, it is quite evident that determining a time lag can be prob- 
lematic. However, we proceed with the usual tools developed for reverberation-mapping 
to see whether a statistically significant detection can still be culled from the data. To 
measure the time lag between the continuum and the broad part of the H/3 emission line, 
we ern ploy the interpolation cross-cor relation function (ICCF) method of iGaskell fc Sparke 
(119861) and iGaskell &: Peterson! (119871 ) and the discrete correlatio n function (DCF) method 



Edelson fc Krolikr ( 1988 ). with the modifications discussed by White fc Peterson (1994) 



of 

in both cases. To qu antify the uncertain ties in the time delay nieasurement, we follo w the 
method described b vlBentz et al. (l20Q6]j) , originally outlined by iPeterson et al.l (119981 ) with 
the modifications of lPeterson et al.l (120021 ). Figure 1 shows the continuum and emission-line 
light curves and cross-correlation results. 



For the light curves of NGC 5548 shown in Figure 1, we find the cross-correlation 
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centroid occurs at Tcent = 6.4^2 3 days and the peak of the cross-correlation function (for 
physical lag times r > 0) occurs at Tpeak = S.Sllg days in the observed fr a me. F ollowing the 



consideration of the relative merits of Tcent versus Tpeak by lPeterson et al.l (120041 ). we will use 
the value of Tcent corrected for time dilation effects by a factor of (1 + z), Tcent = days, 
in the discussion that follows. 

A measure of the quality of the data can be estimated by analysis of the auto-correlation 
function of the continuum in the top-right panel of Figure 1. The ICCF and DCF methods 
give very different values at t = 0, the peak of the ICCF is 1.0, compared with the peak of 
the DCF, ~ 0.5. This discrepancy arises from the fact that we have discarded measurements 
of the DCF where t is exactly equal to (i.e., simultaneous pairs measured from the same 
spectrum) because they are subject to correlated errors, especially when the amplitude of 
variability is so low. The DCF with zero-lag pairs excluded therefore gives a more realistic 
measure of the correlation of nearby data points with each other because, unlike the ICCF, 
it does not include correlation of the noise in the data with itself. With such noisy data 
and low-level signals in the light curves, it is rather remarkable that we are able to obtain a 
statistically significant time lag between the continuum and H/3 light curves. 

To further quantify the significance of this detection, we have carried out the follow- 
ing Monte Carlo simulations. An optical continuum light curve was generated using the 
characteristics of the po wer density spectrum measured for NGC 5548, P{f) oc /"^'^^ 



( ICoUier &: PetersonI l200ll ). The model continuum light curve was then convolved with a 
transfer function to produce a model emission-line light curve. Two simple models for the 
transfer function were employed for illustrative purposes: 



• A thin, spherical shell of radius R = 6.3 light days, which has a constant response from 
T = to T = 2R/c and is zero everywhere else. Such a transfer function results in the 
most smoothing of the continuum light curve shape as it is transferred to the line light 
curve. 

• A thin disk of radius R = 6.3 light days and inclination 0°, which is a delta function at 
T = R/c. This transfer function results in the sharpest transfer of structural features 
from the continuum light curve to the emission-line light curve. 



For ideal data, each of these cases would yield a lag measurement of t = R/c. 

After creating the model light curves, they were sampled in the same pattern as the 
observations presented in this paper. The sampled points were rescaled according to the 
excess variances and mean fractional errors determined for the real data set and presented 
in Table 2. To simulate noise, each data point was randomly altered by a Gaussian deviate 
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within the flux uncertainties appropriate for tlie observing campaign. Tfie simulated set of 
observations was then cross-correlated in the manner described above, and the peak and 
centroid of the cross correlation was recorded. This entire process, beginning with the model 
continuum generation, was repeated 1000 times for each of the transfer functions described 
above. The resulting cross-correlation peak and centroid distributions for the two transfer 
functions are displayed in Figure 2. Casual inspection of the distributions in Figure 2 shows 
that most of the simulation results pile up around the expected lag measurement of 6.3 days, 
which lends much confidence to the measurements made above for our observing campaign. 

About 15% of the simulations result in unphysical lag times, that is, r < 0, indicating 
that the broad line flux responds before the continuum flux varies. Considering only those 
simulations that have physically plausible lag times (r > 0), we find that for the thin shell 
transfer function, 53% of the peak measurements and 56% of the disk measurements fall 
within the range of la uncertainties quoted above for the lag measured in our observing 
campaign. The percentages increase to 69% and 67%, respectively, for the thin disk transfer 
function. The asymmetric shapes of the peak and centroid distributions result in slightly 
lower median lag times than expected, however the medians are well within la of the input 
lag time of 6.3 days: Tpeak = 5.8]'^3;2 and Xcent = 5.5^1'^, based on the simulations with a 
thin shell transfer function; Tpeak = ^-^-l^ and Tcent = ^-^lll'l, based on the simulations with 
a thin disk transfer function. The fact that t he thin disk model seems to mor e accurately 



reproduce the observed data is not surprising: iHorne. Welsh, fc Peterson I ( 119911 ) find strong 



evidence ruling out a spherically symmetric geometry for the H/3 BLR of NGC 5548. The 
most likely geometry is probably somewhere in between these two simple models, however, 
the models are useful in that they approximate the two opposite extremes of possible BLR 
behavior in AGNs in the sense of minimal (face-on disk) versus maximal (spherical shell) 
smoothing of the line light curve shape. In any case, the simulations clearly confirm the 
detected restframe lag of 6.3 days. 



3.2. Line Width Measurement 



The mean and root-mean-square (RMS; i.e. variable) spectra of NGC 5548 were cal- 
culated using the full set of MDM data and are shown in Figure 3. To measure the width 
of the broad emission line, we first subtract the [O III] A4959 narrow line using the [O III] 
A5007 emission line as a t emplate and the standar d scaling factor of 0.340 for the A4959 line 
relative to the A5007 line ( Storey &: Zeippen 200ol ). We also subtract the narrow component 
of the H/3 line usin g the same tena p late a nd the scaling value of 0.110 for the ratio of the 
lines determined by iPeterson et al. J2002h . We then interpolate the continuum underneath 
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the broad H/3 emission line by setting continuum windows on either side of the hne. The 
width of the hne is typically characterized by its full width at half maximum (FWHM) and 
by the second mome nt of the line profile — the line dispersion (crune) — as described by 
Peterson et al.l (120041 ) . and is corrected for the spectral resolution. In this case, since the red 
side of the mean profile can be contaminated by Fe II multiplet 42 emission and artifacts of 
the [O III] removal, we also include in our line measurements the line dispersion based on the 
blue side of the line, assuming the profile is approximately symmetric around the expected 
central wavelength of the line, which we denote crbiue,sym- The uncertaintie s in the line width 



measurements are determined using the method described by lBentz et al.l (l2006bl ). however, 
the formal uncertainties for the cxiine and abiue.sym measurements in the mean spectrum of 
0.6% are probably misleading since the measurement errors are mostly systematic in this 
low state of variability. We therefore adopt an error of 20%, which is more typical of the 
errors in mean profile widths. 

Table 3 lists the line widths measured for NGC 5548 in the mean spectrum, where 
the signal is much stronger, as well as the RMS spectrum. For this particular data set, 
the FWHM of the variable (RMS) broad line is not well characterized, so we omit this 
measurement. Fortunately, all measurements of the line width tabulated in Table 3 are 
consistent with each other even though the RMS spectrum in particular has an extremely low- 
level signal. While the line width measurement tha t is typically preferred for reverberation- 
mapping data sets is (Jiine(RMS) (jCoUin et al.ll2006l ). due to the low-level signal in the RMS 
spectrum, we will choose to use a measurement from the mean spectrum in this case. The 
blending problems between the red side of the H/5 line and the [O III] A4959 line and other 
features cause us to adopt crbiue,sym(niean) as the most credible line width for this data set in 
the following discussion. While we have noted that all the line width measurements for this 
data set are consistent with each other, so the choice of line width does not have a critical 
impact on the following discussion, there are issues that effect certain measurements and not 
others, so a combination of the line width measurements is not appropriate here. 



3.3. Black Hole Mass 

The mass of the central black hole is determined by 

Mb„ ^ (2) 

where r is the time delay of the emission line, AV is the width of the emission line, c is 
the speed of light, and G is the gravitational constant. The factor / takes into account 
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the unknown inclination, geometry, and kinematics of the BLR. lOnken et al.l (120041 ) have 
found that an average value (/) = 5.5 ± 1.7 if the relationship between black hole mass and 
stellar velocity dispersion (Mbh — cr*) for AGNs is normalized to that of quiescent galaxies. 
With this normalization, we find Mbh = ^.Olgy x IO^Mq. Combining this new measurement 
with all previous reverberation-based mass measurements of the black hole in NGC 5548, 
we find a weighted mean of Mbh = 6.83lg;28 ^ lO^M© based on H/3 measurements only, and 
Mbh = 6.541q 25X IO^Mq when including all measurements of broad lines. The reverberation- 
based mass is in good agreement with the mass expected from the Mbh ~ relationship of 
Tremaine et al.l (120021 ) . which predicts Mbh = 9.41^ 2 x lO^Mo based on the stellar velocit y 
dispersion measurement for NGC 5548 of a^, = 183 ± 27 km s~^ (IFerrarese et al.ll200ll ). 
Considering the low-level signals of both the time lag and the line width in this low-luminosity 
data set, it is rather remarkable that our results are consistent with the expectations based 
on previous monitoring programs as well as other methods of estimating the black hole mass 
for NGC 5548. 



4. DISCUSSION 

With the addition of this lowest-luminosity data set to the compilation of NGC 5548 
monitoring data, we examine several relationships between luminosity and other parameters 
where this experiment extends the range that may be explored. 



4.1. Relationship Between Line Width and Time Lag 



Figure 4 shows the width of broad emission lines (crune) measured in NGC 5548 versus 
their time lag relative to the continuum at 5100 A {jcent)- If the BLR kinematics are domi- 
nated by gravity, we expect that there will exist a virial relationship between the time lag and 
the width of an emission line. The top panel shows the relationship for only those measure- 
ments of the H/3 line in NGC 5548, and the bottom panel includes all broad lines available. 
In either case, the slope is basically the same, ab out —0.5 ±0.1, consistent with the virial 
relati onship previously reported for NGC 5548 (jPeterson fc Wandell Il999l : iPeterson et al. 
2004( 1. 
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4.2. Relationship Between Time Lag and Luminosity 

Photoionization models of the BLR are typically parameterized by the shape of the 
ionizing continuum incident on the BLR gas, and an ionization parameter U, defined as 



U 



Q{H) 



(3) 



where Q{H) is the number of ionizing photons emitted from the continuum source per second, 
r is the size of the BLR, and tt-h is the BLR particle density. We would naively expect that 
as the ionizing continuum luminosity varies, each emission line should have the greatest 
response at some particular value ofUnu. This would then lead to the expectation that, for 
any given line, r oc L 



1/2 



In Figure 5 we plot the relationship between the H/? time lag relative to the continuum 
and the luminosity of the AGN at 5100 A. If we include all 13 years of data from the 
AGN Watch campaign and the new low-luminosity data supplied in this work, we find a 
power law slope of 0.73 ± 0.14. However, the data from Year 12 (2000) of the AGN Watch 
program is the most poorly sampled of all the monitoring campaigns for NGC 5548. The 
cross-correlations of the H/3 emission line with the continuu m in the Year 1 2 dat a yields 
very ambiguous results, as can easily be seen in Figure 2 of iPeterson et al.l (120021 ). If we 
remove the Year 12 data point based on its potential unreliability, we find a power law slope 
of 0.66 ± 0.13, within ~ la of the naive expectation of 0.5. Using this relationship and 
the mean luminosity of NGC 5548, we estimate that the lag time for the Year 12 campaign 
sh ould have been r neasu red at ~ 11 days rather than Xcent = 6.6l3;8 days as determined 
by iPeterson et al.l ( l2004l ). With the large error bars due to the ambiguous results of the 
cross- correlations for Year 12, the expected value of 11 days is within la of the measured 
value of 6.6 days. Based on these arguments, we are more inclined to accept the slope of 
0.66 ±0.13. 

One would expect that a measurement of the AGN luminosity in the ultraviolet (UV) 
would give a much better estimate of Lion than in the optical part of the spectrum. In 
Figure 6, we examine the relationship between 83 pairs of optical flux measurements at 
5100 A, Fopt (5100 A), and UV flux measurements at 1350 A, Fuv (1350 A), for NGC 5548. 
The UV and optical flux measurements in each pair were made within one day of each 
other, and no statist ically significant time delay exists between the optical and UV data (see 
Peterson et al.ll2002l for a discussion of the data). The op tical fluxes have be en corrected 
using the host-galaxy starlight contribution measured by iBentz et al.l (l2006a] ). The best 



power law fit to the data is Fopt oc 



0.84±0.05 
UV 



If we combine this with the above r — Lopt fit 
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excluding the data from Year 12, we find r oc L^y Interestingly, this is consistent not 
only with the naive expectation of how an object s hould behave when it varies over time, 
but also with the relationship recently determined by lBentz et al.l (l2006al) for the pop u lation 
of reverberation-mapped AGNs as a whole. On the other hand, ICackett fc Hornd (120061 ) 
find a much shallower slope when modeling the 13 years of AGN Watch data for NGC 5548 
as one continuous light curve with a luminosity-depende nt delay map, rnore in agreement 
with the predictio ns of detailed photoionization models by lKorista fc GoadI (120041 ) . However, 



Cackett &: Hornd used a smaller correction for the host galaxy starlight contribution to the 



flux, which would serve to artificially flatten the relationship they find. 



4.3. Black Hole Mass Estimates 



The most obvious relationship left to examine is whether the black hole mass determi- 
nations are changing with the luminosity of the AGN. We expect that this would not be the 
case, as the black hole mass should not be changing in a measurable way and it certainly 
should not depend on the luminosity. Figure 7 shows the distribution of virial products 
(the black hole mass without the average scaling fact or of 5.5, i.e. , Mrh //) as a function 
of luminosity. This is similar to the recent analysis by ICollin et al.l (120061 ) in that the open 
circles are those based on the measurement of o"biue,sym from the mean spectrum, and the 
filled circles are based on the measurement of aunc from the RMS spectrum. The luminosity 
measurements have been c o rrected for the contribution from host-galaxy starlight using the 
corrections of iBentz et al.l (l2006al ). We have added the lowest lu minosity data point de- 



scribed in this work to the set previously examined by ICollin et al.l . but we use cxbiue.sym from 



both the mean and RMS spectra because of the blending issues discussed earlier in §3.2. 

Inspection of Figure 7 shows that the virial product for NGC 5548 is consistent with 
being independent of luminosity state. This is as expected, since the black hole mass should 
not be varying w ith luminosity. The l argest outlier is again from Year 12 (2000) of the AGN 
Watch program (jPeterson et al.ll2002l ). 



Formal likelihood analysis of the black hole mass measurements gives a width of 0.15 dex, 
which is comparable to the standard deviation of the measurements themselves. From this, 
we may conclude that the intrinsic scatter dominates the errors for repeated reverberation- 
based mass measurements of a single object. 

While the reverberation-mapping experiment appears to be quite repeatable with a 
fairly high precision, we do not have an independent measurement of the systematic errors 
present in reverberation-mapping techniques. In order to estimate the systematic errors. 
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we must compare results with a complementary method of measuring Mbh- Comparison 
of reverberation-mapping results to those derived from the Mbh — cr* relationship shows 
that there is a fac tor of ~ 3 scatter in Mbh measurements from reverberation-mapping 
(jOnken et al.ll2004l ). It is therefore important that inclination and other effects be considered 
in order to reduce the systematic errors in reverberation-mapping results. 



5. CONCLUSIONS 

We present the lowest-luminosity (-L5100 = 4.7 x 10^^ ergs s~^) monitoring data for 
NGC 5548. We measure a rest-frame time lag of 6.3l^2'3 days for changes in the H/3 line flux 
relative to changes in the continuum flux. Combining this new data point with previously 
collected data, we flnd the weighted mean of the black hole mass to be Mbh = 6.54!^q25 x 
IO^Mq. 

With the results of this low-luminosity monitoring campaign, we are able to extend 
the range over which we may examine several relationships that depend on luminosity. We 
conflrm the existence of a virial relationship between the time lag and the line width of the 
broad emission lines in NGC 5548. We flnd that the relationship between luminosity and 
time lag in NGC 5548 is consistent with a virial expectation, and is also consistent with the 
relationship determined for the entire population of AGNs with mass measurements from 
reverberation mapping. 
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Table 1. Continuum and H/3 Fluxes for NGC 5548 
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Table 1 — Continued 



JJJ 


r X [oWv A ) 


rip A40oi 


b ser vat or y 


/OA Kr\ nnn^ 
l^-Z,40U,UUU ) 


ergs s cm A ) 


(10"^"^ ergs cni~^) 


uocie 


3465.949 


6.58 ± 0.19 


2.51 ± 0.08 


M 


3466.070 


7.04 ± 0.09 




H 


3466.922 


7.33 ± 0.21 


2.52 ± 0.08 


M 


3467.973 


6.52 ± 0.19 


2.50 ± 0.08 


M 


3468.957 


6.45 ± 0.19 


2.55 ± 0.08 


M 


3469.922 


6.45 ± 0.19 


2.42 ± 0.07 


M 


3470.012 


6.49 ± 0.09 




H 


3470.926 


6.53 ± 0.19 


2.48 ± 0.07 


M 


3471.516 


6.73 ± 0.32 


2.31 ± 0.04 


C 


3471.922 


6.34 ± 0.18 


2.46 ± 0.07 


M 



^The Julian Date listed is the midpoint of the observation 

^Observatory Codes: C = CrAO 2.6-m Shajn Telescope + Nasmith Spectro- 
graph; H = Haleakala Observatories 2.0-m MAGNUM Telescope + MIP; M = 
MDM Observatory 1.3-m McGraw-Hill Telescope + CCDS. 



Table 2. Light Curve Statistics 



Sampling 
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1.39 ±0.14 
1.42 ±0.06 



^Fluxes are in the same units as those used in Table 1. 
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Table 3. H/5 Line Width Measurements 
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Mean 


RMS 


Width 
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2662 ± 532 
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6396 ± 167 





'^The preferred measurement of the 
width of H/9 for the new data presented 
in this paper. 
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Fig. 1. — The left panels show the light curves for the continuum region at 5100 A and for 
the H/9 A4861 line. Measurements taken within a 0.5 day bin have been averaged together. 
The flux is in units of 10~^^ ergs s^^ cm^^ A^^ for the continuum and units of 10~^^ ergs 
s^^ cm~^ for H/5. The shaded area of the continuum light curve shows the contribution to 
the flux from the host galaxy. The right panels show the result of cross-correlating each 
light curve with the continuum light curve; the top right panel is therefore the continuum 
auto-correlation function. The sohd hne shows the ICCF method, while the points show the 
DCF method, as described in the text. 
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Fig. 2. — Cross-correlation peak distributions (solid lines) and centroid distributions (dashed 
lines) for the Monte Carlo simulations described in the text. The top panel is for a thin shell 
transfer function, and the bottom panel is for a thin disk transfer function. 
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Fig. 3.— Mean and RMS of the MDM spectra in the observed frame of NGC 5548. The RMS 
spectrum still shows contributions from the [O III] A4959 and A5007 lines due to imperfect 
subtraction of all the nightly line profiles. The spike at 4825 A appears to be due to a bad 
column in the detector. The dashed lines in both windows show the wavelength limits and 
the dotted hues show the continuum levels used in calculating the width of the broad H/3 
line. The flux density for both spectra is in units of 10"^^ ergs s~^ cm~^ A~^. 
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Fig. 4. — Rest-frame emission line widths versus time lags for NGC 5548. The top box 
shows the relationship for all measurements of the H/3 line, while the bottom box shows the 
relationship for measurements of all lines. The filled circles sh o w the measurement of H/? 
from this work, while the open circles are from iPeterson et al.l (120041 ). The solid lines are 
the best fits to the data, and the dotted lines are fits with a forced slope of -0.5, i.e. a virial 
relationship. 
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Fig. 5. — Rest-frame H/ 3 time lag versus the luminosity at 5100 A. The open circles are 
from iBentz et al.l (j2006al ) after the luminosity measurement has been corrected for the host 
galaxy starlight contribution, and the filled circle is from this work, also corrected for host 
galaxy starlight. The dashed line has a slope of 0.73 ±0.14 and shows the fit to the data with 
all the points included. The solid line shows the fit excluding the data from Year 12, with 
a slope of 0.66 ± 0.13, which we take to be the best current fit to the data for NGC 5548. 
The dotted line has a forced slope of 0.5, as one would expect from naive predictions of the 
behavior of the gas in the BLR to changes in the ionizing flux. 
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Fig. 6. — Optical flux at 5100 A, after correction for host galaxy starlight, versus the ultra- 
violet flux at 1350 A. Each pair of optical and UV measurements was made within one day. 
The solid line shows the best fit, which has a power law slope of 0.84 ± 0.05. 
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Fig. 7. — The virial product determined with o"iinc calculated from the mean spectrum (open 
circles) and the RMS spectrum (filled circles) versus the luminosity at 5100 A after correction 
for the contribution from host galaxy starlight. 



